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Abstract 



We construct a natural invariant measure concentrated on the set of square-free num- 
bers, and invariant under the shift. We prove that the corresponding dynamical system is 
isomorphic to a translation on a compact, Abelian group. This implies that this system is 
not weakly mixing and has zero measure-theoretical entropy. 
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Introduction and Notations 

Let V be the set of prime numbers. By p (with or without indices) we will always denote 
an element of V. A positive integer n is square-free ii \ n for every p. Denote the set of all 
square-free numbers by Q (for quadratfrei) . The indicator of the set Q is the function n ^ ^"^{n), 
where fj, is the Mobius function: 



The functions fi and fi'^ are of great importance in Number Theory because of their connection 
with the Riemann zeta function. For example. 



Furthermore, the estimate |X]n<7v /^('^) | ~ O^^N^^'^^'^) as — t- oo is equivalent to the Riemann 
Hypothesis. P. Sarnak [S] has recently addressed a number of statistical and ergodic properties 
of the sequences (/i(?T,))„ and (/^^(n))^. 
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if n 



if n is not square- free; 

if n is the product of k distinct primes. 



1; 
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0.1 Notations 

We shall use the standard notation e{x) = e^'^'^. For every integer n denote by u{n) the number 
of its distinct prime factors. For example, w(l) = and uj{2 ■ 3) = u{2^^ ■ 3^) = uj{7 ■ 23) = 2. We 
shall also use the notations 

V{n) = {p : p\n}, V2{n) = {p : p'^\n}. 

Notice that if n G Q, then = u}{n), V2{n) = 0, and V2{n'^) = V{n). For every finite set 

A<ZV, define 

[A] = \{p. 

p&A 

In particular [0] = 1. Notice that if A, B are disjoint, then [AUB] = [A] [B] and [Ar}B\ = 1. 



1 Formulation of the results 

The goal of this paper is to describe a dynamical system 'naturally' associated to Q and study 
its statistical and ergodic properties. 



1.1 Correlation functions 

The first step is the construction of correlation functions for Q. Choose r integers < /ci < 
k2 < . . . < K and consider the set 

QNih, k2, . . . , K) = {n < N : /i^(n) = fi^{n + ki) = . . . = fi^{n + kr) = 1}. 

The ratio 

n 1 I \ \QN{ki,k2,...,kr)\ , . 

EN{ki,k2,...,kr) := — (1) 

is the frequency of square- free integers n < N for which n+ki,n+k2, . . . , n+kr are also square- free. 
It also gives the expectation (hence the notation E) of the product iJ?{n)ij?{n + ki) ■ ■ ■ ji^{n + kr) 
with respect to the uniform measure on {1, 2, . . . , A^}. Notice, by taking r = 1 and ki = 0, that 
Qjv(O) is simply the set of all square-free numbers not greater than A^. It is well known that 

6 

lim Ejv(O) = ^ ^ 0.6079271018 (2) 



We include the proof of ^ and some of its generalizations in Section [2| see Theorems 2.1 2.3 



The study of E,i\f{ki, kr) as N ^ oo is also classical, see L. Mirsky R.R. Hall P, K.M. 
Tsang [H]. It is known that the limits 

Cr+i{ki, . . . , kr) = lim En {ki, kr) (3) 
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exist. This fact holds true because the hmit in the rhs of (|3j) is expressed in terms of averaging 
over arithmetic progressions. We shall refer to c^+i as the (r + l)-st correlation function for Q. 
Various formulae for Cr+i{ki, . . . ,kr) are known (see Section El). We prove that 



c.^^ih,...,K)= E E E (-iF-'"' E Fy^- 

0</'<«"<r fj,'2{di,j„) = 1 mo,mi,...,mr>0 ViCV,0<l<r ^U/=0 ' 'J 

df, i„\ki, - kii, \Vi\=mi 

[Vi, nVi„] = di,^i„ 

To the best of our knowledge, the formula Q is new and, although complicated, plays a role in 
the spectral analysis of the correlation functions. Let, for example, r = 1. For every d E Q define 

mo,mi>0 -p^^-p-^c^-p L 

|Pol = "10, |Pi| = mi 
[Po r\Vi\=d 

Explicit formulae for ad are given in Section |4j Then 

C2{k)= E ^'^ (6) 

/.2(d) = 1 

and the corresponding spectral measure v on (i.e. satisfying C2{k) = v{k) by Bochner theorem) 
is pure point, given as sum of 5-functions at the points e{l/d?), where d E Q. More precisely, 

^= E ^d^5e(l/d^), (7) 

^l^^(d)=l 1=0 



where the convergence of the series is guaranteed by Lemma 3.2 The spectrum (i.e. the support 
of z/) is the group 

A = |e(^^^ : 0<l<d^-l, ii\d) = fi\gcd{l,d^j) = l| C §^ 

Notice that every element of A is represented uniquely. Moreover, every rational number of the 
form ^ such that d is square-free, < Z < rf^ — 1, and gcd(Z, d"^) is also square free can be written 
as 

— = — ++— (8) 

rf2 pI p2^' 

for some h, . . . , Im, where {pi, . . . ,Pm} = V{d). This representation (fsl) is unique if one imposes 



the restriction < Ij < pf^ — 1, 1 < j < m. In other words, the group A is isomorphic to the 
direct sum ^^Z/p^Z (where only finitely many coordinates are non-zero). Therefore, A is the 
Pontryagin dual of the direct product group 



G = JJ Z/p% (9) 



p 
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which is an Abehan compact group (endowed with the product topology). In other words, G = A. 
Each element g G G is identified with a sequence {gp2)p^-p indexed by V, where gp2 e Z/p'^Z: 

g = (5'4,5'9,fi'25,5'49, • • •)• 

Given h G G, denote by Th : G — ?■ G the translation Th(g) = g+h. Let B be the natural a-algebra 
on G, and let us put the uniform measure on each Z/p^Z. The corresponding product measure P 
on B is invariant under translations, and therefore it is the Haar measure. The ergodic properties 
of translations on compact Abelian groups were studied for the first time by J. von Neumann. 
He showed [12] that two ergodic such translations with the same spectrum, are isomorphic as 
measure-preserving dynamical systems. This is true in general for ergodic transformations with 
pure point spectrum and it plays an important role in our analysis. Later, P.R. Halmos and 
J. von Neumann P] proved that every ergodic dynamical system with pure point spectrum is 
isomorphic to a translation on a compact Abelian group. This implies, for example, that every 
ergodic dynamical system with pure point spectrum is isomorphic to its inverse. For an historical 
survey on the isomorphism problem see 0. 



1.2 A Natural Dynamical System 

Consider the space X of all bi-infinite sequences x = {x{n), — oo < n < oo} where each x{n) 
takes value either or 1. Denote by B the natural cr-algebra of subsets of X, and introduce the 
probability measure n defined on B as follows: For every r > and every — oo < ko < ki < . . . < 

kr < oo 

n {x G X : xik^) = x{ki) = . . . = x{kr) = 1} = Cr+i{ki — ko, k2 — k^, . . . , kr — ko), (10) 

where c^+i is the (r + l)-st correlation function (|4]). It is clear that (10) determines the measure 
n uniquely. We call 11 the natural measure corresponding to the set of square-free numbers. 

If T is the shift on X, i.e. Tx = x', x'{n) = x{n + 1), then it follows immediately from (10) 
that n is invariant under T. We can now formulate the main result of this paper: 

Main Theorem. 

(i) The dynamical system {X, B, 11, T) is ergodic and has pure point spectrum given by A. 

(ii) (X, B, n, T) is isomorphic to (G, B, P, T^), where u = (1, 1,1,...). 

P. Sarnak [9] formulates the result that (G,B,P, Tu) is a factor of {X,B,Il,T). His methods 
also allow to show that the factor map is in fact an isomorphism. Our approach is rather different 
and is based on a spectral analysis of the dynamical system {X,B,Il,T). The statement in the 
following corollary can be also found in [9]. 

Corollary. The dynamical system (X, B, H, T) is non weakly-mixing, and its measure-theoretic 
entropy is zero. 
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It is worthwhile to remark that the main focus of [9] are the topological dynamical systems 
M = (0^(„),T), S = (0^2(„),T) given by the shifts on the orbit closures of {^{n))n and (yU^(n))n, 
respectively. The topological entropy of S is positive, equal to ^log2, and R. Krishna in his 
senior thesis |3] identifies a measure of maximal entropy for S. 

Our paper is organized as follows. Section |2] includes the classical computation of the density 
of square-free numbers and its generalization to square-free numbers avoiding finite sets of prime 
factors (the proof is given in Appendix |A]) . The latter will be used for the computation of certain 
relevant constants. Section [s] contains the proof of formula (|6]) for the second correlation function 
(see Theorems 3.1 an d |3.3[ ), and formula (|4]) for arbitrary correlation functions (see Theorems 3.7 



and 3.10). In Section we compare our formulae for the correlation functions with the ones by 
L. Mirsky and R.R. Hall. We also present several useful lemmata (some of which are proven in 
Appendix |B|. The spectral analysis for the dynamical system (X, i3, 11, T) is carried out in Section 
[5] and yields the first part of our Main Theorem. The analysis of the spectrum for (G, B, P, Tu) 
is done in Section |6} and the second part of our Main theorem follows from it, by means of a 
theorem by J. von Neumann 
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2 The density of Q and some of its subsets 

Recall that E7v(0) = ^ |{n < iV, n G The following theorem is very classical. 
Theorem 2.1. 

lim Ejv(O) = (11) 

Proof. We can write /i^ as the indicator of the set of square-free numbers by imposing the condition 
that its argument avoids all arithmetic progressions modulo p^: 

^?{n) = \{{l-xA^)) (12) 

p 

In the above expression Xp2(^) is the indicator of the arithmetic progression {p^l : / G Z}. Let 
us open the brackets in ( [l2| : 

/i2(n) = I - ^Xp-^in) + ^ Xpi{n)Xpl{n) - Xpl{n)Xpl{n)Xpl{n) + • • • • 

P Pl<P2 Pl<P2<Pi 
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We can write 

n<N p ^ pi<p2 pi<P2<P3 

_ 1 _ 6 

P^J C(2) vr2 

Here and below denotes a remainder that tends to zero as — ?■ oo. □ 

The statement of Theorem 1 can actually be refined as follows: 
Theorem 2.2. 

En{0) = — + 0{N-^/^) asN^oo. 



In other words, e^v in the proof of Theorem 2.1 satisfies the estimate Ien] = 0{N^^^'^). This 
result is also very classical, and is a special case of Theorem 2^ below. Let us fix a finite set 
S <ZV and define the set 

2:^(0) = {n<N : fi^n) = 1, peS ^p]n} (13) 

of all square-free numbers not bigger than N and not divisible by any of the primes p E S. For 
example, Q\f (0) is the set of odd square-free numbers not bigger than A^. Notice that when S 
is empty we get the full set of square-free numbers, i.e. 2^(0) = Qn{0). In analogy with ([T]), let 
us define 



N 



We have the following 

Theorem 2.3. For every finite S G V we have 
where 

and the constant C{S) implied by the Og-notation can be taken as 

pes ^ \pes / pes 

The proof of Theorem 2.3| is presented in Appendix|Aj it implies the existence of the asymptotic 
densities 
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For example, the density of the set of odd square-free numbers is (i.e. odd and even square- 
free numbers are in 2:1 proportion). Analogously, by choosing S = {p}, we see that the set of 
square-free numbers not divisible by p is "p times as large" (in the sense of density) as the set of 
those divisible by p. If, for instance, we choose S = {2, 3} we obtain a({2, 3}) = 1/2, and we see 
that 50% of the square-free numbers is not divisible by either 2 or 3. 



3 The Formulae for the Correlation Functions 

We first address the formula for the second correlation function C2- Recall that Ejv(fc) 
jj: \ QN(yk)\. We prove the following 

Theorem 3.1. Let k be square-free. The limit 

C2{k) = lim Kiy{k) 

N~^oo 



exists and 



Mk)= E (-1)— E py^- (15) 



mo,mi>0 Po,^! C V 

{Vol = mo, IVil = mi 
VonVi=0 



Notice that the quantity in the rhs of (15) does not depend on k. Moreover, Vo^iVi = <^=^ 
[Po n Pi] = 1 and [VoflVif = [Vq U Vif. This motivates the definition (g of a^. Furthermore, 
in ^ we can add the constraint mi,m2 > uj{d) in the outer sum, as this is necessary in order for 
[Vq n Vi] = d to hold for sets Vo,Vi of cardinality mi,m2 respectively. We prove the following 
formula for 0"^: 



Lemma 3.2. 



^^ = in(i-;l)- (16) 



p\d ^ ^ 

Proof. Recall that, since d is square-free, 1^(^)1 = uj{d). By setting m = mi — uj{d) and M 
mi + m2 — 2uj{d) in ^ we obtain 



<rm<r A.r V / -n/ ^ -n . 



d^lV']^ 

0<m<M ^ ' -pi (^-p ^ -p((i) i 



\V'\ = M 

2 



Ly y (zar^in 1 

^2 r-p/12 ^2 111-' p2 

M>0 V'cV^Vid) ^ J p\d ^ ^ 

\V'\ = M 



□ 
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In particular, Lemma 3^ shows that cr^ is positive and bounded away from zero and infinity. 
More precisely 

6 



< (Ti < (Td < 



where cti = Hp (l - ^) ~ 0.3226340989. We can rewrite 

p\d ^ 

The general formula for the second correlation functions is given by the following 
Theorem 3.3. Let k be an arbitrary integer. The limit 

C2{k) = lim E,N{k) 

exists and 



C2 



o- Cl- 



in) 



P-Hd) = 1 



Proof of Theorems 3.1 and 3^. Let k be square-free. 



n<N n<N ^ p[ v'x<v'i 



Xp'^{n)Xp'i{n)Xp'i{n) + • • • ' 1 " XI ^p'i"(^ + ^) + 



Pi 



+ V + ^)^P2" + ^) - Xpf + A;)Xpf + ^)Xpf {n + k) + ... 



p'i<p'i 



After opening the brackets in (18) we can re-write it as a sum of expressions of the following 
form: 



p"i <p",<--<Pi" 

'•1 *2 '"m\ 



(19) 
(20) 



of prime numbers with cardinality mo and m\ respectively. Observe that the product in (19) can 
be written as 



where n< mi,m2 > 0, and Vq = {p[^, 1 < / < mo}, Pi = {p'l^, 1 < I < mi} C V are subsets 

,nd mi respectively. Obser^ 

Xp'.2 (^)Xp'2 (n)--- Xp'2 (n) = X[PoP (^) • 

*1 ^2 *"io 



In the same way 



Xp//2(n + k)Xp>>^{n + k) ■ ■ ■ Xp'H^ {n + k) = Xir^rin + k). 
The product X[Va]'^ij^) " X[Pi]^{.'n + /c) is 1 if and only if there are two integers Iq, h such that 



n 

n + k 



Then ( [2lp2| imply that 



k = [V,]% - [Vo]%. 



(21) 
(22) 

(23) 



The products [Vq] and [Vi] are both square-free. Denote by d their greatest common divisor 
(gcd). Then (23) has solutions IqJi only if k is divisible by (P. We can reformulate the last 
statement differently. Choose a square-free d and consider pairs of sets Vo,Vi C V such that 
gcd([Po]5 [Pi]) = d (or, equivalently, [Vq fl Vi] = d). By dividing by both sides of (23) we have 



k 
'd? 



(Po) 



d^ 



-h- 



d^ 



By definition of gcd, it is clear that {ViY/d'^ and {VoY/d'^ are coprime integers. Thus, one can 
find two integers ao, ai such that 



rf2 



c/2 



-ao 



and all possible values of Iq^Ii in (|21[ 22|) have the form 



k 



From (21) 



k 



c/2 



(i2 



k 



ai 



rf2 



[n]'^«o+lcm([Po],[Pi])'/, 



(24) 



where 1cm is denotes the least common multiple. Equation (24) means that n belongs to an 
arithmetic progression with step lcm(['Po]5 ["^i])^ = [Po UPi]^. Given d, 1711,1712 and Vi,V2, after 
summation over n and dividing by N, we get l/lcm([Po]5 [Pi]Y + ^n- 

If k is square-free, then the only square-free d that divides k is d = 1. This means that Vo 
and Vi are disjoint (i.e. [Vq fl Vi] = 1 and [Vq U Pi] = ['Po]['Pi]), and otherwise arbitrary, finite 
subsets of P. In this case 



mo,mi>0 'Po.'Pi C P 

|Pol = "10, l^il = mi 

-Po n Pi = 



(25) 
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This concludes the proof of Theorem 3.1 For arbitrary k we have 



c,{k) = lim Mk) ^ J2 H 

Ai2(d)=l,d2|fcmo,mi>0 -pg^-p^cV 

IVol = mo, \Vi\ = mi 
[Vo nVi] = d 

/l'^(d) = l,d2\k 



[^0UPi]2 



and this gives Theorem 3.3 



□ 



Note that Theorem 3.3 includes Theorems 2.1 and 3.1 as special cases. In particular from ^ 
and (17) we get 

' (26) 



E 

/.2(d) = l 



The same result follows from the explicit formula ( 16 ) for a^- 



Remark 3.4. Theorems 3.1 and 3.3 show that the value of C2{k) depends on the arithmetic 
properties of k. This fact is certainly very unusual from the point of view of Probability Theory 
and Statistical Mechanics. If k is square free, then the function C2{k) takes the constant value 
(Ti. Analogously, C2{k) is constant along any subsequence of numbers k sharing the same set of 
divisors that are the square of a square-free number. If we define T>{k) := {d : = 1, d'^\k}, 

then = T^{k') =^ C2{k) = C2{k'). The opposite implication follows from the formula ( [l6| . 

Observe that every set 'D{k) is of the form 



V{k) = \\{p:V'<Z V2{k) \ 



(27) 



where V2{k) = {p : p'^\k}. This means that \V{k)\ = 2l^2(fc)| ^nd V{k) = V{k') V2{k) = 

V2{k'). The set of k such that 7^2 (^) = is the set of square-free numbers, and we know that it 
has positive density (equal to G/vr^, given by d2])). In general, we have the following 



Proposition 3.5 (Density of the level sets of C2). Fix a square-free number d. Then the density 
of those k's such that C2{k) = C2{d'^) exists and is given by 



^{d') := hm h{k<N: C2{k) = C2{d')}\ = \ 

p\d ^ 



{2i 



Proof. If V{d) = {pi,P2, ■ ■ ■ ,Pm} then k satisfies C2{k) = C2{d'^) if and only if it is of the form 
k = p1^ ■ ■ ■p'^q, where /i^(g) = 1, > 2, and pj f q for every j = 1, . . . ,m. Fix Oi, . . . , a^, > 2. 
Then 



\{k<N, k=p\^ ■■■p'^^q, /^'(g) = 1, p.fgfor J = l,...,m}| 



1 



-E 



{pi,...,Pm} 



(0) 
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and, by Theorem 2^, the hmit as — )■ oo is 

6 1 



TT 



n =TT - 



Now, by summing over all Oj > 2, we obtain 



6 



TT — V — = -TT - = -TT — 



and the proposition is proven. 
Remark 3.6. We can check that 



□ 



m>0 -p' c P pGP 



1 + 



p2 _ 1 



6 



-1 



/.2(d) = l 

Here we present the values of ^((i^) for square-free numbers d < 17. The sum of the corre- 



sponding densities is ~ 97.6% and one can check that X](i<42- f (rf ) > 99% 



1 


2 


3 


5 


6 


7 


10 


11 


13 


14 


15 


17 


6 


2 


3 


1 


1 


1 


1 


1 


1 


1 


1 


1 


^2 


7r2 


47r2 


47r2 


47r2 


87r2 


127r2 


207r2 


287r2 


247r2 


327r2 


487r2 



d 



One can also compute the limit 



1 ^ / 
lim — > Cofn) = I 

n=l ^ 



0.3695753612 



by considering the series ^^2(^)=i d{d'^)c2{d'^) and using Proposition 



3j| 



and Lemma 



3.2 



We shall 
structure 



retrieve this fact from the more general result of Lemma [4. 5[ Figure [IJ summarizes 
of the second correlation function. 

Let us now consider higher order correlation functions. Recall ([T]) and set ko = 0. We have 
the following 

Theorem 3.7. Let ki,...,kr be such that all the differences kii — kin, < I' < I" < r are 
square-free. Then the limit 

Cr+i{ki, A;2, . . . , kr) = lim EN{ki, /c2, • • • , K) 



exists and 



C-r+l {ki, k2, ■ . . , kr 



E 



E 



mo,mi,...,mr>0 



Vi CV, < I < r 
Vi, nVi„ = 0forl' 7^ I" 



(29) 
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•••••• - C2(6^) = 0-1+0-2 + 0-3+0-6 

- • • " C2(2^5^) = 0-, + 0-2 + 0-5 + o-|o 
C2(2^) = 0-1+0-2 

• - C2(3^5") = o-| + 0-3 + 0-5 + 0-15 

. . . ... ... ... ... ... . . - C2(32) = 0-1 + 0-3 

. •-c2(52) = o-i +0-5 

— — . — — - ..... — . - — — .V ... — . .. - — — — .. - C2(l) = (T, 

50 100 150 200 250 

k 

Figure 1: The second correlation function C2{k) and its level sets. 



Remark 3.8. Notice that the rhs of (29) depends neither on ki,...,kr nor on the values of 
their differences as long as they all are square-free. This shows that c^+i does not decay as its 
arguments grow. In other words, the set of square-free numbers has arbitrarily long correlations 
of size 

Remark 3.9. It is not enough to check that the consecutive differences ki—ko, ^2—^1, • • • , kr—kr-i 
are square-free in order for all differences to be square-free. For example, if (/ci, ^2, ^3, ^4) = 
(1, 6, 7, 10), all consecutive differences are square free but 2'^\ki — ^2 and 3^|A;4 — ki. 

The general case of correlation functions is addressed by the following 

Theorem 3.10. For general ki,k2, . . . k^ the limit 



Cr+i{ki, k2,..., K) = lim E^iki, /c2, ■ ■ ■ , K) 

N^oo 



exists and 



Cr^,{k„...,kr)= E E E (-iF-"^' E FT-^- ^30) 

0</'<«"<r /i2{d,, ;„) = 1 mo,mi,...,mr>0 Vi, < I < r lU/=0 ' 'J 

dli„\ki,-ki„ \Vi\=mi 

[Vi, nVi„] = di,^i„ 
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Proof of Theorem 3/l_ and 3.10 . We have 



n<N 



n<N 1=0 



^ -^X(^pmj2{n + ki) + X^pW)2{n + ki)X(^pW^2{n + ki) + 

\ 

X^pWyiin + fci)X(p(0)2(n + ki)X(^pm-j2{n + h) + . . . 



J 



where /cq = 0. After opening the brackets we have a sum of expressions of the form 



■1) 



mQ+m\+ . . .+inr 



n 5Z X(p{0)2("- + fc/)---X(p(o )2(n + /i;/). 



'=0 p« 



Since for each /, the primes pf^ ,pf^ , . . . ,pf^^ are all distinct, the product YYjlLiX(^pW)2 is the 



indicator of the arithmetic progression with the step YlTliiPi^'')'^^ 



n 



+ h = Uipfjfsi, < / < 



(31) 



All the sets {pf\ 1 < J < rrii} are pairwise disjoint. In fact, if we assume that the sets P'^''^ = 
{Pi. , 1 < J < "^i'} and ' = {p\. , 1 < j < fni/'} have non-trivial intersection, then ki' — kin 
should be divisible by fl P*-'"-*]^, contrary to our assumption that kii — kin is square-free. As 
in_the proof of Theorems |3.1 and 3^, for every pair I', I" we can describe all solutions si', si" of 



(31). Find integers s), , s,„ such that 



a')x2 fo) 



1 = 11(P 



ii J ^l' ' 



i=i 



Then 



Si' = {kv - ki»)sf^ + sY\_{p. 



ii J ' 



J = l 
m,, 

Si" = {ki' - kv')sf) + s^{p. 



(n\2 
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From (31) we have 



n 



(V)\2 



j=i j=i j=i 

for arbitrary t. In the same way 



(32) 



n 



('')^2TTUnA2 

i ' 

j=l j = l j = l 



(33) 



with the same t. It is easy to check that the rhs of (32) and (33) coincide. Furthermore, (32) and 
(33) show that n belongs to an arithmetic progression with step ni=o YY^iiPij)"^ ■ Therefore 

lim K]\r{ki, k2, ■ ■ ■ ,kr) = Cr{ki, k2 — ki, . . . , kr — fcr-i) = 

= y (_i)ELo-. y — I — 

mo,mi,...,mr>0 "pC) C "P, | = m; l'^l=0 ' J 

0<l<r 



This expression is a direct generalization of (25) and concludes the proof of Theorem 3.7 



Consider now the general case of arbitrary ki, 1 < I < r. Formula (31) means that n belongs 
to an arithmetic progression with the step UTliipl^f- Fix 

a square- free number c?;/ ;// such that 

ki' — kill is divisible by d^iin- Then the products YYjliipf^Y^ ^ — ^^■■■^^^ should satisfy the 
condition 



vj=i j=i 

This gives the possibility to write the final expression for Cj.+i(/ci, . . . ,kr): 



Cr+l{kl, . . . , kr) = X] XI XI *^ 

0<l'<l"<r /^2(d;,_,„) = 1 mo,mi,...,mr>0 

dl in\kv - kin 



-Pj, < « < r 
\Vl\=mi 

[Vi,nVi„] = dii^i„ 



and Theorem 3.10 is thus proven. 



□ 



Remark 3.11. Notice that Cr+i{ki, . . . ,kr) might be zero if r > 3. For example, 04(1, 2, 3) = 
since there is no n such that n, n+1, n+2, n+3 are all square- free. All cases when Cr+i{ki, . . . , kr) = 
correspond to constraints modulo for some prime p. This fact is clearly reflected by the general 
formula for Cr+i(/ci, . . . , K), which is discussed in Section 111 
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4 Other formulae for the correlation functions 

L. Mirsky [1] proved that 

Ap \ki, . . . , kr) 



Cr+l{ki, k2,...,kr) = Y\il 

V \ 



(34) 



where Ap'^^\ki, ^2, • • • , ^r) = |{0, A;i(mod p^), fc2(mod p^), . . . , kr{mod p^)}|. Notice that 

l<A'^+'\k^,k2,...,kr)<r 

for finitely many p and Ap^^\ki, kr) = r + 1 for infinitely many p. For r = 1, we have 



1 

2, otherwise. 



This gives, for instance, 



- n (i - ,4) n (i - ,4 



(35) 



We can check that the formula (17) for C2{k) coincides with (35). Indeed, since V2{k) = {p : p'^\k} 
and V{k) = {UpeV'P- ^' ^ ^2(A;)}, we have 



C2{k) 



nO-^)n 



^ ^ defCit) p\d ^ 



pz 



pz 



ev{k) p\d 

n 



pZ 



P&V2{d) 



(36) 



In particular, if /c = 0, then 2^(0) = Q and 7^2(0) = V and thus (36) allows us to retrieve (26): 



E 

M2(d)=l 



(37) 



The case of Ap~^^\ki, . . . ,kr) = r + 1 corresponds to the case when 0,ki, . . . , kr are distinct 
modulo p"^ for every prime p. This means that the differences fc/' — ki^r are not divisible by p"^ for 
every prime p. In other words, the differences ki/ — kin are all square-free. This case is addressed 
by Theorem 3.7 Let us check that the formula (29) for Cr+iiki, k2, . . . , kr) corresponds to (34) 
with A^~^^\ki, . . . ,kr) = r + 1 for every p: 
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m>0 -p' -p I- J p 

|P| = m 



Analogously, one can check that the formula (30) is equivalent to the general formula (34) with 



no restrictions on ki, . . . ,kr. R.R. Hall P proved the following useful 
Lemma 4.1. For every < ki < k2 < ■ ■ ■ < kr we have 



Cr+i{ki,...,kr) = ^---^g{so)---g{sr) ^ e(%ki + ... + ^kA , (38) 



]>1 Sr>l < < - 1 

M2(gcd{t„s2)) = l 
< j < r 

+ ^ + ■ . ■ + |i e ; 



where ^ 

(lis) = — TT 



^(^) = ^M^)n^- (39) 



Moreover, the series in (3S) converges absolutely 



4.1 Averages of the correlation functions 

We formulate here three interesting lemmata, based on the above formulae. They concern the 
average of the second correlation functions along arithmetic progressions and will be used later. 
The proofs are given in Appendix [B| 

Lemma 4.2. Let d be square-free. Then 

1<N ^ ' peVid) ^ 

Lemma 4.3. Let d be square-free and letl <t < d'^ — 1, gcd((i^, t) = g >1, where g is square-free. 
Then 

i<N ^ ' pev{d) ' 

Lemma 4.4. Let d be square-free, and let 1 < t < d"^ — 1, gcd(t, d"^) = g > 1. Then 
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The following two lemmata deal with exponential sums involving the second and the third 
correlation functions. Recall the function g from (39). 

Lemma 4.5. Let A = e(^) G A. Then the limit 



1 ^ 

2)2(A)= hm ^$^A"C2 



n] 



n=l 



exists and 2)2 (A) = g{dY . 

Proof. We can write n = d'^l + 1 for some I > and < t < rf^ — 1 and set 



(40) 



where 



Lemma 4.2 



gives 



C2{d^m + 1). 



m<{N~t)/d? 



lim (A) „ 



1 / 6 A ^ TT- p2 



For t 7^ 0, the value of \im.N^oo I^nW given by Lemmata 
where g = gcd(t, d"^). More explicitly, 



4.3 



4.4 



(41) 

It depends only on V2{g), 



lim I^'\X) = e 



lt\ 1 / 6 
d?) '(PK^ 



p\d p^\g 



(42) 



Let us introduce the notation 



n 



p2|gcd(t,d2) 



p2-2' 



Notice that ri((i^) = r^i^tid'^) and therefore the limit lim7v-!.oo -?^7v(A) is real. Using (41) and (42) 
we can write 



lim /at (A) 

N^oa 



6 



p\d ^ \ p\d ^ t=l 



27Tlt 
~d?' 



rt{d^ 



(43) 
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and, if u{d) = 17^(^)1 = r, then 



cos 



cos 



pi,P2\d t < - 1 
Vlvl\i 

- E E 

Pl,P2,P3\d. t<d?-l 



2TTlt 
27111 



1 - 



E E 

Pi\d t<d'^-l 
Pl\t 

pl-r 



cos 



27llt 



pI-2 



+ 



pJ-1 

pI-2J\' Pi -2 



cos 



rf2 



1 - 



Pi-2J\ Pi -2 



pi,...,Pr^i\d t < d'^ - 1 

p?---p?-llt 



COS 



Pl-l 

p\-2 



P^-i - 2 



+ ...+ 



(44) 



Recall the assumption that gcd(/, d^) is square-free, and notice that for every pi, . . . ,PmM, m < r, 



E 

i < d2 _ 1 



COS 



-1. 



Now (44) yields 



-1 -1 



E ) ^'C^) - (-^' E ,-5^ E ,5 _ _ , 

*=1 ^ ^ V PiM^^ Pi,P2\d^^ ^2 



^ -1 -1 -1 

^ pI-2pI-2pI-2 

Pl,P2,P3\d 



■'y-' E ^ 



pi,...,p,._i|d 



- 2 pi_, - 2 



n - - Mr n,-f^) - n,-^ - n 

^ ^ / ^ ) p\d'^ p\d 



p^-i 

p2 -2' 



and (43) becomes 



lim /Ar(A) 



(I) i+n^-i =(|) \[j^--o(d) 

^ ^ p\d^ \ p\d^ I ^ ^ p\d > 



Remark 4.6. Since jp^ <p^ — l<p^ for every p, then 



-<n— 

p\d ^ 



< 



□ 



(45) 



Since d is square-free, if we want to give an upper bound for u:{d) in terms of as d — )■ cxd, it is 
enough to consider the case when d is the product of the first r prime numbers: d = piP2 ■ ■ - Pr- 
In this case u{d) = r. It is known that logc? = rlogr(l + o(l)) as r — )■ oo. This means that 
in general u{d)loguj{d) < (1 + £:i)log(i for every ei > 0, provided that d ^ 1. This implies 
uj{d) < ^(°og'^^i+ei) ; where W denotes the Lambert function, i.e. the solution of the equation 



X = W{x)e^^^\ It is known that W{x) ~ logo; as x — t- oo. Therefore uj{d) < 
every 62 > and thus 



logfji+^i 
log(logdl+=i)l-=2 



for 



uj{d) 



Oeidn 



(46) 



for every e > as d ^ 00. In other words, formula3(41 42) give 



as (i — )■ 00, 



for every t = 0, 1, . . . , c?^ — 1. However, the cancellations coming from the different exponential 

2-irit I I 

factors e~d^ in lN{d) are responsible for the higher order of smallness shown in (43): 



lim lN{l,d') = 0, 



as (i — )• 00. 



Lemma 4.7. Let \i,\2 G A, Ai = e^^j , A2 = <2(^^j, and X = A1A2 = e(^) G A. Then the 
2-fold limit 



2)3(Ai,A2) = hm -i- K'K'c3in,,n2) 

Ni 00 I\il\2 



N2 ^ 00 



ni=l n2=l 



exzsts anc? 2)3 (Ai, A2) = g{di)g{d2)g{d) . 



Proof. Using Lemma 4.1 we can write 



2)3(Ai,A2) = .hm ^YYXTXTYT.J29i^Msi)9is2) 



Afi -> 00 N1N2 
N2 —>■ 00 

E 



ni=l n2=l 



3>1 Si>l S2>1 



^1 ^2 



i = 0,1,2 

Sp -^1 ^^2 



(47) 



(4J 



Let us bring the limit and the sums over ni,n2 in (48) inside the sum over to,ti,t2- For fixed 
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So, Si, S2 we have 

E 



< < , 



1 



1 

lim — > e 

ni = l 



iVi^oo iVi 



+ 



dl 



rii 



J = 0,1, 2 



N2 

lim } e 

Af2-s>oo 

n2=l 



/2 



n2 



The two sums over ni , n2 can be written as 



S 



2 rf? 



r 1 - +^ 

1, 



if 0^1 



where = e(j^ + and j = 1,2. Thus, as — )■ oo, only the indices tQ,ti,t2 such that 
^1 = ^2 = 1 give a non-zero contribution to (48). This condition means 



s? rf? 



1, ifA, ^1; 
0, if \j = 1; 



for j = 1, 2. However, because of the conditions < < s^ — 1 and /i^(gcd(tj, s|)) = 1, the index 



^{d]-l,), ifA,^l 
0, ifAj = l 



will be considered only when s,- = dj. The value of to is given consequently by 



|. ifAi^l = A2 



if Ai = 1 7^ A2 
0, if Ai = 1 = A2. 



In all cases this means ^ = ^, and the condition /i^(gcd(to, Sq)) = 1 implies that sq = d and 
to = In other words ,(48) becomes 2)3(Ai, A2) = g{di)g{d2)g{d), and the lemma is proven. □ 

A2 = ^ ^ satisfy gcd{di,d2) = 1, then 



Remark 4.8. Notice that if Ai = e(^^j 

2)3(Ai,A2) = 



6 



1) 



u}(d\d2) 



n 

p|did2 



1 



(p2 _ 1)2- 
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The product HpidlP^ ~ 1) appears in several formula above. Concerning this product, we 
have the following basic 



Lemma 4.9. Let d he square-free. Then 



- 1) = 1(1 < / < rf' - 1, ^^\gcA{l,dP)) = 1}\ 

p\d 



Proof. By standard inclusion-exclusion we can write the rhs of (49) 



Pi\<i Pl,P2\d Pi,P2,P3\d 



-I). 



p\d ^ ^ ^ p\d 



4.2 Spectral analysis of C2 



We can rewrite (17) 



C2{k)= J2 ^dik), where Ka{k) :-- 

/.2(d) = l 



0"^; if (P\k] 
0, otherwise. 



(49) 



□ 



(50) 



The function is constant (equal to aa) along the arithmetic progression {Id'^ : / G Z} and 
elsewhere. This function is the Fourier transform of a measure on the circle S^, given by a sum 
of 5-functions at the points e{l/d'^), Z = 0, 1, . . . , — 1, with equal weights 0^1 (P ■ A corollary of 
this fact is the formula ([7| for the spectral measure z/ on 



5 The spectrum of the shift operator T 



Recall the definition of the dynamical system (X, i3,n,T) given in Section 1.2 Denote by U 
the operator on the Hilbert space "H = i^^(X, i3, 11) given by 



(t//)(x) = /(Tx). 



(51) 



Since T is measure-preserving, the operator V is unitary. The goal of this section is to prove the 
following 

Theorem 5.1. The spectrum of the operator U is given by A. 

Let us show that A is contained in the spectrum of U. This fact is given by the following 
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Theorem 5.2. Let X = e(^) G A. Then there exists a function 9x E H such that 

{Uex){x) = Xex{x) (52) 

for U-almost every x E X. 

Proof. Let /o(x) = a;(0) and let U\ be the unitary operator on l-i defined by 

{Uxh){x) = X~^h{Tx). 
By the von Neumann Ergodic Theorem, the following limit exists in l-L: 

N N N 

lim -yf/-/o(x)= lim -V A-"/o(T"x) = lim - V A'^xH =: ^a(x). (53) 



n=l n=l ri=l 

The function Q\ is f/A-invariant, i.e. {U\Q\)[x) = 6x{x) for H-almost every x E X. This implies 

□ 



that X-^OxiTx) = 9xix), i.e. (52). 



Denote by x{s) the function X {0, 1} given by the projection of x G X onto its s-th 
coordinate. Introduce the subspace H CTi, 



H = |^asx(s)|, 



i.e. the closure of the set of all complex linear combinations of the x(s)'s. H is invariant under 



U, and by (53), all the eigenfunctions Ox belong to H. Let us remark that, since the operator U 



is unitary, the eigenfunctions 9x are orthogonal to one another for different A. Let us write 

^^s) = Y,{x{s),ex)ex. 

AeA 



Recall (39). We have the following 



Theorem 5.3. Let X = e(^) E A. Then for every s e7L we have 

{xis),ex) = ygidy. 



Proof. Let us use (53) and write 



{x{s), Ox) = lim ( x(s), ir7 V A 



x[n) 



1 ^ 

lim — \^X'^{x(s),x(n)). 



(54) 



(55) 



n=l / n=l 

Notice that {x{s),x{n)) = C2{n — s), where C2 is the second correlation function given by Theorem 



3.3 Equation (55) becomes 



TV N 

{xis),ex)=hm ^A«C2(n-.) = AMim ^ A"c2(n) = A^ 2)2(A). (56) 



n=l 



n=l 



The needed statement follows now from Lemma [4.51 



□ 
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Theorem 5^ immediately implies the following 
Corollary 5.4. All eigenf unctions Ox are non-zero. 



Remark 5.5. The formula (53) can be written for arbitrary measure-preserving map, but in 
most cases (e.g. automorphisms with continuous spectrum) it gives zero. Theorem 5.3 shows 
that in our case it is non-zero. 

We can also compute the L^-norm of each eigenfunction explicitly. 

Theorem 5.6. Let A = e(^) G A. Then 



\0y 



W)\. 



(57) 



Proof. This is a straightforward application of Theorem |5.3 

1 ^ 



\\o> 



x[n) 



1 ^ 



n=l 



n=l 



n 

p\d 



(p2 _ 1)2- 



□ 



Proposition 5.7. The set of eigenf unctions {Ox}xeA is a basis for H. 



Proof. Since the eigenfunctions are orthogonal it is enough to show that they span the space of 
all linear combinations of the a;(s)'s. We know that each atom {A} of the spectral measure u 
(associated to the second correlation function via Bochner's theorem) corresponds to Ox in the 
space H generated by linear forms, and these form a set of generators for H. □ 

Let us define the normalized eigenfunctions: for A G A set 

Ox 



0, 



\0, 



so that {Ox}x&A is an orthonormal basis for H. Let us write 



AeA 



Since {^a}a is an orthonormal basis for H, then by Lemma 



4.9 



and Theorem 

2 



5.6 



n 1+ 



deQ p\a 

6 



r'cv,\v'\<oopev 
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The same argument allows us to estimate the size of the error term in the following approximation 
of x(s): for D > 1 let 



A = e{l/d'^) e A 
d< D 



Arguing as in Remark 4.6, we have 



\x{s) - xois) 



A = e{l/d^) G A 
d> D 



x(s), 9) 



2 6 



J2W)\ = 0,{D-'^n (58) 



d>D 



for every e > 0. 

Let us consider the product of two eigenfunctions Ox^ and 9\^. We have the following 
Theorem 5.8. Let Ai = e(^^j , A2 = e(^^j E A. Then 

dXidX'z = f ^A, 

where A = e(^) = A1A2 and e = e(Ai, A2) = fi{di) ix{d2) iJ,{d) = ±1. 
Proof. It is enough to show that for every s G Z we have 

'9jx,,x{s)) = (ex,x{s) 



Using the definition (53) we can write 



Ni N2 

e,,e,, = hm K'^'K^'^^Mxin,) 

A^i -> 00 I\il\2 — — 



N2 00 



and thus 



n=l 71=2 



Ni N2 



{9,,e,„x{s)) = hm -^J^Y. {x{m)xin2),xis)) . 



A^i -s> 00 N1N2 
N2 ^ 00 



n=l n=2 



Notice that {x{ni)x{n2), x{s)) = 03(^1 — s,n2 — s). Therefore, by Lemma |47f 

A^i N2 



{ex,9x2,x{s)) 



lim 

A^i 00 Ai A2 

N2 00 



71=1 n=2 

Ni N2 

A- hm Tr^EE^r"^^r^^3K,n2) 



TVi 00 A"! A'2 
A''2 -> 00 



n=l n=2 



A-^2)3(Ar\A2') = X'' 9id,)gid2)gid). 



(59) 
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On the other hand, by Theorem 5.3 



{ex,x{s)) = e-'^X)2{X) = X-'g{df 

Therefore 

g{di)g{d2)g{d) \g{d)\ 



\g{di)\\g{d2)\ g{dy 

□ 



By associativity of muhiphcation, e(Ai, A2)e(AiA2, A3) = e(A2, A3)e(Ai, A2A3). Theorem 5.8 



can 

be apphed iteratively. It allows us to write all polynomial expressions in the eigenfunctions as 
linear expressions, and this is a very important fact. 

We want to show that the set of eigenfunctions {^aIasa is a basis for the whole space H. We 
shall need the notion of unitary rings introduced by V.A. Rokhlin (see [7]). 

Definition 5.9. A complex Hilbert space is called a unitary ring if and only if, for certain 
pairs of elements, a product is defined satisfying: 

(I) if fg is defined, then fg = gf, 

(II) if fg, {fg)h and gh are defined, then {fg)h = f{gh), 

(III) if fh and gh are defined and a, /3 G C, then (a/ + (3g)h = a{fh) + P{gh), 

(IV) there exists e G such that ef = f for every / G i^, 
(V) if fng are defined and fn ^ f, fng — ^ h, then fg = h. 

(VI) The set 971 = {/ G : fg is defined for all g E Sj} is dense in S); moreover if fg is defined, 
then there exist /„ G 9Jt such that fn^f and fng fg, 

(VII) for every f E S), there exists f E such that {fg, h) = {g, fh) for all /, (? G VJl. 

An important result by Rokhlin is that every unitary ring can be written as = L'^{M, Ai,m), 
where {M,Ai,m) is a Lebesgue space (see, e.g., V.A. RokhliE[^ [6J). In our case we have the unitary 
ring "H = L^(X, i3, n) and the subspace H which is a sub-ring because of Theorem 5.8 In this 
representation a subring C corresponds to a cr-subalgebra of A/f, i.e. = L^{M,Af, m\j^). 
Therefore if is a subspace of "H, which is a Hilbert space corresponding to some a-subalgebra J-" 
of B. Let us show that 

Proposition 5.10. Up to sets of measure zero, = B. In other words, H = Ti. 



^The notions of Lebesgue space used here allows point with positive measure, contrarty to the classical case 
discussed by Ya.G. Sinai [TD] in the context of K-systems 
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Proof. Let us use the technique of measurable partitions by Rokhhn (see [8]). According to it J-' 
corresponds to some measurable partition ^ of X. \{ C. B, then there exists a bounded, non- 
negative function h{x) and a subset A G J-" such that E(/i|C^) > a for almost all & A and some 
positive a. As any measurable function h can be approximated arbitrarily well in L°°(X, J-", n|jr) 
sense by a function /i' which is a polinomial in the x(s)'s. Using (58) we can approximate h' in the 
measure sense by a finite polynomial in the eigenfunctions 9\. However, every such polynomial 
belongs to our Hilbert space L'^{X, J-", n|jr) and it is measurable with respect to J-". Therefore the 
conditional expectation of h' with respect to ^ is arbitrary close to /i', but such a function cannot 
approximate h in measure. This shows that H = 1-1. □ 



Propositions 5.10 and |5.7| immediately give the following 



Corollary 5.11. The set of eigenfunctions {^aIasa o basis in the space T-L. 



This fact, together with Theorem 5.2 and Corollary 5.4 yields Theorem 5.1 and the following 



Theorem 5.12. The dynamical system (X, i3, n,T) is ergodic. 

Proof. By constructing the eigenfunction 9i (Theorem 5.2) and showing that it is non-zero 
(Corollary 5.4), we proved that the eigenspace corresponding to the eigenvalue 1 is at least 



one-dimensional. On the other hand, by Theorem |5.1[ its dimension cannot be bigger than one. 
This implies that the only invariant functions are constants H-almost everywhere, and hence we 
have ergodicity. □ 



Theorems 5.1 and 5.12 give part (i) or our Main Theorem. 



Remark 5.13. One could also derive Corollary |5.11| in a different way and without using the 
theory of unitary rings and measurable partitions by Rokhlin. The derivation, although explicit, 
is rather long. In fact, one can show that for every —oo<Si<...<Sr<oo the product 
x{si) ■ ■ ■ x{sr) belongs to the span of {^aIasa- For example, for r = 2, by Theorems 5.3 and 5.8 



x{si)x{s2) 



Ul6A / \A2GA 

Ai,A2eA 



E 

AeA 



M Yl KK'9{di)g{d2) 



Ai, A2 e A 

A1A2 = A 



and one can prove that 



Y W,)g{d2)\ = OM 



-2+e\ 



Ai, A2 e A 
A1A2 = A 
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for every e > 0, where A = e(^). This imphes that 



E 

AeA 



J2 K'K'9idi)9{d2) 



Ai, A2 e A 

A1A2 = A 



is finite. 



6 Spectral analysis for (G, B,P, Tu 

Recall the group G defined in ^ . Let us consider the space EI = 
operator U on EI defined by 

(U/)(g) = /(g + (1,1,1,...)). 



and the unitary 



>.. It is immediate to see that the 



Theorem 6.1. The spectrum o/U is given by A. 

Proof. Consider the projection TTp2 : G Z/p^Z, 7rp2(g) = 
function ^e(i/p2)(g) = e(7rp2(g)/p^) is an eigenfunction for U with eigenvalue ^(^^^- By taking 

powers one can get any eigenfunction ^e{t/p2) with any eigenvalue ^(^^j < t < — 1. By 

multiplying different such eigenfunctions (with different p), one can obtain eigenfunctions with 
an arbitrary eigenvalue A G A. Since A is the character group of G and is a translation in G, 
then there are no other eigenvalues. □ 

To conclude the proof of part (ii) of our Main Theorem we need the following 



Theorem 6.2 (J. von Neumann, |I2])- Two ergodic measure- preserving transformations with 
pure point spectrum are isomorphic if and only if they have the same spectrum. 



Theorems 6.1 and 6.2 imply that (X, i3, n,T) and (G,B,P, Tu) are isomorphic as measure- 
preserving dynamical systems. This concludes the proof of our Main Theorem. 



A The Proof of Theorem 2.3 



This Appendix is dedicated to the proof of Theorem |2.3 It is based on the following identity 

(a * b)(n) 



b{n) 



vra=l / \n=l 

where a * 6 is the Dirichlet convolution of a and b: 



E 

n=l 



n" 



(a*6)(n) = 5^a(rf)6(f) 

d\n 



(60) 



(61) 
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We shall be considering only the case of s = 2 and bounded sequences a{n) and b{n), therefore 
there will be no question about convergence of the above series. We shall also use the classical 
identity 

^/i(d) = 0, forn>l. (62) 

d\n 

First, let us consider the case of square-free numbers not divisible by a single prime p, i.e. 
S = {p}. In this case, we shall prove Theorem 2.3 by means of three lemmata, and then we shall 



explain how to generalize this approach to general finite sets S. 

be the indicator of the integers not divisible by p, i.e. 



Let Wp{n) 



Wp{n) 



0, if p\n; 



otherwise. 



Lemma A.l. 



fi^{n)wp{n) 



J2 Mwp{d)wp{^) 

d: d2|n 



(63) 



Proof. If p\n, then p\d or p\^ (possibly both) for every divisor d of n. Thus Wp{n)wp[^) = and 
the sum in the rhs of (63) is (and obviously equals the Ihs). If p f n, then no divisor d of n 



will be divisible by p and Wp{d) = Wp(^)=l. The sum in (63) then becomes "^d^in f^i^) ■ ^ 
is square-free, then c?^ = 1 is the only perfect square that divides n, and the sum equals 1 (and 



clearly agrees with the Ihs of (63)). If n is not square-free let us write it as ni?T,2 where rii and 
77,2 are defined as follows. For every prime p let us define I = l{p,n) = max{j : p'ln}] then set 



77,1 



'[lp,2^iP and n2 = 11^:21;^^ ' IIp:2hP ' • ^^^^^ square free, if d' 



Lp: 2\l ' 



V:2\l 



then d\n2. This 

means that the sum in the rhs of (63) becomes ^^^i^^ t^i.^) and equals by (62) (thus matching 
the Ihs). This concludes the proof of ^he Lemma. □ 



Lemma A. 2. 



n=l 



p^ 



Proof. The formula follows from the trivial computation 



E 

n=l 



Wp{n) 



^ 77,2 ^ 



n=l 



C(2). 



□ 



Let us denote by {Si{n)}n>i the sequence equal to 1 if 77 = 1 and otherwise. Then 
Lemma A. 3. 

{^Wp) *Wp = Si. 
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Proof. For n = 1 the statement is obvious since c/ = 1 is the only divisor of n and we have 
/i(l)wp(l)^ = 1 = 5i{l). Let n > 1. Then {(f^Wp) * Wp) (n) = X]cz|n • 
discuss the cases when p\n and p f separately, and argue as in the proof of Lemma A.l[ In the 
first case we have that Wp(n)wp(^) = and the sum is 0. In the second case Wp{d) = Wp(^)=l 
and the sum becomes J2d\nf^i^)' ^^^^ t)y (62). In other words, we have shown that, for n > 1, 
we have((/iWp) * Wp) (n) = = Si{n), and this concludes the proof of the Lemma. □ 



Corollary A. 4. 



E 

n=l 



IJ,{n)wp{n) 



p 



p2-iC(2)' 



□ 



Proof. This corollary is a straightforward application of Lemma A. 2 and the formulae (60, 61) 
with a = fiWp, b = Wp, and (from Lemma A. 3) a*b = 6i. 

We can now give the 

Proof of Theorem 2.3 when S = {p}. Notice that |Qj^"'^(0)| = X]n<Af By Lemma 
we can write 

'^fi\n)wp{n) = '^'^fi{d)wp{d)wp{^) . 



A.l 



(64) 



n<N 



n<N (P\n 



Now we want to exchange the two sums. Let us fix c? < yN . For every n of the form n = md"^ 
we have Wp(^) = Wp{m)wp{d). Let vlt^N) be the number of integers of the form md^ where 



< ^ and p\m. Then 



\Q^S\^)\ = E v'f\N)^^{d)wp{d). 

d<\/N 



We can estimate the number vi][^ [N) as follows. Let [^J = t (mod p), t G {0, 1, 2, ... ,p — 1}. 
Then 



where 



q\'\d,N) 



< 



p — 1 
p 



p 

N 
d^ 



(p-l)+t 



This gives us 



where |g2^''^(^)l < C'{{p})VN. Now, Corollary 

p 1 



Q^\0)\ = N 



N 

P — 1 
p 



t l- 



Np-1 
p 

p — 1 
p 



+ qi'\d,N), 



< 2 



p 



p 



-■■ cm). 



d<VN 



A.4 



\Q]^\n)\ 



p+lC(2) 



yields 

iV + gW(iV) + gW(Ar), 
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where 



p — I 

p — ' ct^ p 



dx 



p-l N 



P VN -1 



< C"i{p})VN 



for every > 4, where C"{{p}) = 2^—. This concludes the proof of the theorem, with a{{p}) 
^ and Ci{p}) = C'iip}) + C"i{p}f. 



□ 



Let us now discuss how to adapt the above proof of Theorem 2.3 for the case of a general 
finite set of primes S C V. The sequence Wp has to be replaced by the indicator of the integers 
divisible by none of the primes in S, i.e. 



ws{n) 



0, if p\n for some p E S; 

1, otherwise. 



Lemma A.l is still valid if we replace Wp by ws- 



Lemma A. 5. 



d: d2|„ 



(65) 



Lemma A. 2 is replaced by an analogous statement given by inclusion-exclusion: 
Lemma A. 6. 



E 

n=l 



ws{n 



where 



<s) = n 



a(5)C(2), 

p2 



pes 

Proof. If 5 = {pi,P2, ■ ■ ■ ,Pk}, then inclusion-exclusion gives 

k 



a{S) 



=1 ^' l<ji<i2<A: 



E 

l<ii<i2<«3<fc 



{PhPi2y {Pi,Pi2Piz)'' 



+ . . . + 



(-1)^ 



(piP2 ■■■pky 



i=i ^ 



□ 



Lemma |A.3| also holds for ws '■ 
Lemma A. 7. 



ifiws) *ws = Si. 
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Finally, Corollary A.4 is replaced by 
Corollary A. 8. 

^ fi{n)ws{n) 

n=l 

We are now ready to give the 



a(5)C(2) 



Proof of Theorem 2.3 for general S = {pi, . . . ,pk} C V. Lemma [A. 5 gives 



n<N (P\n 



(66) 



Let us fix (i < \/N . For every n of the form n = md^ we have ws{^ = ws{m)ws{d). Let ri^{N) 
be the number of integers of the form m(P where m < ^ and p f m for every p E S. Then 

d<VN 

The set of numbers not divisible by any p E S, has density given by 



i=l l<h<^2<k 



1 -1 
— + ... + — 



1=1 l<li<i2<fC l<il<l2<?3<fe 



1=1 



pes 



The estimate of rj'^iN) comes from the following observation. If 



^2 



= t (mod [S]) for t G {0, 1, 2, . . . , [S] - 1}, 



then 



where 



[S] 



pes 



pes ^ 



p — 1 
p 



< 



< 



pes ^ \pes / pes 



This gives 



\Q%m = Ni[ 



p — 1 
p 



pe5 d<^/N 
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where |gf (A^)| < C'{S)Vn. Now Corollary IZsl yields 



iQ^(o)i=n 



p 



pes 



p+lC(2) 



N + q^{N) + ql{N), 



where 



pes 



d>VN 



and C"{S) = 2 Ylpes ~ — ^'^^ N > A. This concludes the proof of the general case of the theorem, 
with aiS) = U,es ^^I'and 0(6) = C'iS) + C"{S). □ 



B The Proofs of Lemmata I4.2M4.4 

Proof of Lemma \4.^ Let us write / as follows: 

per p'ev 



(67) 



where V C V{d), a{p) > 2 for every p eV,V' CV \ V{d), \V\ < oo, b{p') > 2 for every p' e V, 
q is square-free and p \ q for every p eVU V. It is clear that every I > 1 can be written uniquely 
as in (67). And the condition / < can be rewritten using the notation in (13) as 



Furthermore, notice that C2{(P'l) can depends only on d and V: 



c2{dH) = n 

peP{d)uP' 



p^ 



n 



p^V{d)VJV' 



1 - 



P^ 



n 



peV(d)vjv 



p^-l 



Now we can write 



^E^^c^o-. n fci E E n'- 

■ peV(d)^ VQV(d)V'CV\V(d)p'eV' ^ 



KN 



12 



12 



- YlpeviPY^^^ ■Ylp'ev'ipT'''^ ^/(n,e^(p-)''('^'-n,'e."(p')''('''')^^^ 



a(p) > 2, p e P 

Kp') > 2, p' e V 
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Now we can use (14) while taking the hmit as — oo, and sum over all a{p), b{p') > 2 as in the 
proof of Proposition 3.5 Notice that the sets V{d) and V are disjoint. We obtain 

i<N pev(d)^ Pcv{d)V'cv-^v{d)p'eV' ^ pePuP' 



- TT -^T\^—=i-\ TT ^ 

T,2 11 p2_lll 2 11 2_1' 



and the lemma is thus proved. □ 

Proof of Lemma \4.3[ Let us first consider the case g = 1- Numbers of the form Ilp'ep' P'''^^ '* ' 
where V <Z V \ V{d), \V'\ < 00, > 2 for p' G P', q is square-free and p f g for every 

p EV{d) UV can be represented as 

Yl .q = dH + h (68) 

for some 1 < < — 1, where gcd(/i, d^) = 1. Since there are '■p{d'^) such /I's (here (^9 denotes 
Euler's totient function) and the various /I's appear with the same frequency, then 

lim 4 y\ C2{dH + t) = ^ lim 4 "E E ^^^^H + /i) (69) 

1<N ^ 1<N gcd{h,d2)=i 



Notice that the condition / < becomes 

'id2N+h)/{Up,^^,p'Hp') 



and 



Now 



cAsi + = n (1 - i) n (1 - 1) = n fc^- 



^E E o,{dH + h) = a, n 



«<Af gcd(/i,d2)=i V'CV^Vid) p'eV ^ 

d'^N + h pi^d)yjv' 

b(p')>2,p'eP' \\.\.p'(^v'P ) 
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and by taking the limit as — )■ cxd we get 

,3 



n ^ E n ^;rr^=^i;^ n ^ n (^+^ 

6 -p-r p^{p^ -2) -p-rp2_l / pi{pi _ 2) 



Let us apply the fact that ip is multiplicative and that ipijP') = p{p — 1). We obtain 



Now (69), (70) and (71) yield the desired result. Let us now consider the case when gcd(t,(P) = 
g = p. In this case (P = p'^df and t = pti, where di is square- free, p \ di, and p\ ti. We can write 

-l + -=pdll + t,= l\p"^^">-qu 

where V C V ^ '^{d), \V'\ < oo, qi is square-free and p f g for every p G V{d) UV. The condition 
/ < reads now as 

^1 ^ -^ipdiN+t,)/{n,,^r'P"'^'''^r 

Since, by assumption, \ d'^l + t, then we have 

c,{dH+t)=c,{pdi+u)= n (i-i)n(i-|)=-- nfcl 

Now, since gcd{ti,pd\) = 1, we can use (69): 

^ E + ^) = ^ E + ^i) = i™. ^ E E + h), 

(72) 



Af^oo ^ N-^oo N ^ '^(pd'^) N^oo N 

1<N 1<N 1^ 1<N gcd{hupdl)=l 



and we can write 

1<N gcd(hi,pdl)=l V'CP-.V(d) p'eV ^ 

Pdf pdlN + hi 'pi^d)vjv' 

b{p')>2,p'(^r' ^^P '^^ ^ ^ ^ 
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Notice that V{d) = V{di) U {p}. By taking the hmit as — )■ cxo we obtain 

i^^^E E C2{pdii + h) = 

1<N gcd{hi,pdl)=l 
J2 



E n n ^ n ^ 

6 p' 

<7i 



1 11 „ I 1 11 „/2 _ 



7r2p+l AJ. „+l ^ J.J.„/2_2 



n ^ n (i+ 



6 — 2 T-r p^(p2 — 2) T-r p2 — 1 

6y p\p'-2) yr p\p^-2) 

Let us use the fact that ^{j)d\) = (p{j>)(p{d\) = {p — l)(p{dl) to obtain the formula 

1 



^ n (74) 



'f(P4) P - 1 p(p 



Now we can combine (72), (73), and (74) to conclude the proof of the Lemma. The case of a 
general square-free g is treated analogously. □ 

Proof of Lemma \4.4\ The case when g is square- free (i.e. V2{g) = 0) is already included in 
Lemma 4.3 Thus, it is enough to consider the case when V2{g) 7^ 0. Let, for simplicity, 
gcd(t, d"^) = g = p^ (i.e. u{g) = 1), the case of u{g) > 1 being analogous. We have that = p^d\ 
and t = pHi, where di is square-free and p f di. In particular gcd(ti, df) = 1. We can write 

where a > 0, V C V \ 'P{d), \V'\ < 00, qi is square-free and p \ qi for every p G V{d) U V. The 
condition I < N can be written as 

and 

peP'u{p} ^ ^ p^P'u{p} ^ ^ peP'u{p} ^ 
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Notice that by V and {p} are disjoint by construction. Using (69) we see that 



4E^2(rf'/ + t)= lim Y,C2{p\dll + h)) 



N^OO N N^OO ■ 

KN KN 



^iL-^E E -'(mi^i^d)- (75) 



We have 



and by taking the hmit as — t- oo we get 



P - 1 V- TT -1 6 p yr 1 

p2_2 2^ 11 „/2 _ 2p- 1 11 r)+ 1 11 

^ V'cv^r{d) p'ev ^ ^ p&r{d) ^ p'ev ^ 



6 p+1 



n^^ E n 



^^-^ ' P&V{d)^ V'cV\Vid)p'&V' ^ 

6 p+1 T-r T-r / 1 



6 P+1 TT TT A 

^V2nfn2_2) p + 1 V 



"7r2p(p2_2) 11 p + 1 11 V p2_2 

^ Q TT - 2 -|-|- p2 - 1 ^ / 6 y p2 - 2) 

7r2p2_l 11 p+lp2_lll 2 1^2) p2_l 11 (p+l)(p2_lV ^ ^ 

We use again the fact that 

1 1 



and combining (75), (76) and (77), we obtain the Lemma. □ 
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